Dynamical trapping and relaxation of scalar gravitational fields 



m 
o 

(N 
C 

(N 



cr 



> 
On 

(N 
O 

o 



x 



C. H.-T. Wang, 1 ' 2 '* A. 0. Hodson, 1 A. St J. Murphy, 3 T. B. Davies, 1 R. Bingham, 2 ' 4 and J. T. Mendoncca 5 ' 2 

1 SUPA Department of Physics, University of Aberdeen, King's College, Aberdeen AB24 SUE, UK 
2 STFC Rutherford Appleton Laboratory, Chilton, Didcot, Oxfordshire 0X11 OQX, UK 
3 SUPA School of Physics and Astronomy, University of Edinburgh, Edinburgh, EH9 3JZ, UK 

y 'SUPA Department of Physics, University of Strathclyde, Glasgow G4 ONG, UK 
'' 'GOLF '/Centro de F'isica de Plasmas, Instituto Superior Tecnico, 1049-001 Lisboa, Portugal 

We present a framework for nonlinearly coupled scalar-tensor theory of gravity to address both 
inflation and core-collapse supernova problems. The unified approach is based on a novel dynamical 
trapping and relaxation of scalar gravity in highly energetic regimes. The new model provides a 
viable alternative mechanism of inflation free from various issues known to affect previous proposals. 
Furthermore, it could be related to observable violent astronomical events, by specifically releasing 
a significant amount of additional gravitational energy during core-collapse supernovae. A recent 
experiment at CERN relevant for testing this new model is then briefly outlined. 
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Relativistic theories of gravity with a scalar field <fi, 
that actively rescales the metric g a b used by Einstein's 
general relativity (GR), have been a subject of fundamen- 
tal significance in theoretical and experimental physics, 
astronomy and cosmology [1]. Such scalar-tensor (ST) 
theories of gravity arise universally from unified gauge 
theories and are important for approaches to quantum 
gravity involving a conformal symmetry [2]. However, in 
view of the stringent experimental constraints [3] within 
the solar system on the linearly coupled ST theory due 
to Brans and Dicke [4, 5], recent works have focused on 
nonlinear theories under extreme gravity. These include 
inflation in the early Universe [6, 7] and strong field ef- 
fects in neutron stars, pulsars and core-collapse super- 
novae (SNe) [8-14]. 

As is well known, the observed large scale structure of 
the Universe, including its flatness and horizon, suggests 
an inflation shortly following the Big Bang. The original 
"old" proposal [15] involves the trapping of a hypothet- 
ical scalar inflaton field in a false vacuum associated by 
a prescribed potential. As the Universe supercools this 
metastable vacuum energy functions as a cosmological 
constant to drive the exponential expansion of cosmic 
scale, until the end of the first-order phase transition to 
the true vacuum state through bubble nucleation. Un- 
fortunately, inflation does not end naturally in this sce- 
nario and its nucleation also leads to unacceptably large 
subsequent fluctuations in density and cosmic microwave 
background (CMB) [17]. 

While a "graceful exit" from inflation could be 
achieved through an alternative "new" slow-rolling scalar 
model [18] without invoking bubble nucleation, at the ex- 
pense of a new fine tuning of the inflaton potential, such 
an option is considered to be ad hoc. Likewise, other 
inflation models such as chaotic inflation suffer from var- 
ious fine tuning drawbacks [19]. 

Remarkably, if ST theory is used to evolve the cosmic 
scale factor a(t) during inflation, then its resulting power- 



law growth, as opposed to the exponential growth in the 
new inflation, would find a graceful exit through the first- 
order phase transition. The first attempt of such an "ex- 
tended inflation" [19] using the linear theory requires a 
Brans-Dicke coupling outside the observational bounds 
[16, 20]. Fortunately, this flaw can be surmounted by 
using general ST theories in "hyperextended inflation" 
[6] , which retains the power-law growth of the scale fac- 
tor a(t) and the associated graceful exit. Its continuous 
development can be found in the recent review [7] and 
references therein. 

Here we raise a new avenue of inflation using general 
ST theories, that not only assimilates the strengths of 
previous models while avoiding their drawbacks but also 
offers a number of attractive features. We shall assume 
the scalar gravitational field to play the role of the in- 
flaton field and postulate no further fields violating the 
strong energy condition. As with the new inflation, a 
potential is required to allow the scalar to roll downhill 
during inflation. The slowness of rolling and the initial 
uphill positioning of the scalar need not be fine-tuned. 
In contrast, it is "controlled" through a novel dynami- 
cal trapping mechanism as a result of a generic nonlinear 
coupling between the scalar and matter. The effective 
vacuum is initially seized in a state with a large cos- 
mological constant and then relaxes to one with a small 
cosmological constant in a dynamically controlled man- 
ner without bubble nucleation. In a sense the model re- 
news new inflation without fine tuning using ST theory. 
Nonetheless it recovers GR in the low energy domain of 
the current epoch. In addition, the dynamical trapping 
mechanism could be reactivated in certain observable en- 
ergetic astronomical events such as core-collapse SNe. 

The new model is constructed from the general La- 
grangian for ST theory [5, 21]: 



g 1/2 Kd*x + £ + L 



16ttG 



2 



in the Einstein frame, where a, b — 0, 1, 2, 3, x° = ct, 1Z 
is the scalar curvature of the "Einstein metric" g a b with 
signature (-, +, +, +), 
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is the Lagrangian for the scalar gravitational field 
with a potential V((j>), and L is the action for matter 
fields coupled to gravity through the "physical metric" 
g a b = r2(0) 2 5 , a&- Similarly, in what follows, we shall re- 
fer to a barred quantity as a "physical" quantity, which 
is defined with respect to g a b as opposed to g a b- The 
conformal mapping uses a coupling function 0(0), with 
a coupling strength between the scalar and matter pro- 
portional to a(4>) — 91nQ(0)/50 The function In O(^) 
depends linearly on in the original ST theory by Brans 
and Dicke [4, 5], but becomes nonlinear in a general ST 
theory. Strictly speaking, the gravitational constant G 
here is defined with respect to g a b- However, anticipating 
9ab <Jab in the late Universe, we shall simply evaluate 
G as usual, and use it to define the Planck length l p and 
time t p in the normal way. 

The field equations for g a b and follow from varying 
the total Lagrangian as follows 
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where G a b is the Einstein tensor and □ is the 
Laplace-Beltrami operator with respect to g a b, T ab = 
2g~ 1 / 2 5L/5g a b is the effective Einstein stress tensor for 
0, T ah = 2g~ 1 / 2 5L/5g a b is the Einstein stress tensor for 
matter, with the trace T = g a bT ab . They are related to 
the physical stress tensor for matter and its trace by 



Tab — O 2 T a b, T — ft T 



(3) 



via the conformal mapping between g a b and g a b- 

We shall now apply the above framework to cosmol- 
ogy with a view to scalar gravity induced inflation using 
parameters relevant for energetic astrophysical events in 
the present epoch. Motivated by the recent interest in 
the quadratically coupled scalar gravity in neutron stars 
and core-collapse SNe [8-14] , we shall consider the second 
order coupling specified by: 

0(0) = e 5^ 2 i.e. a(0) = /30 

for some constant /3. We adopt the scalar potential 

V{4>) = £0 CT + A /2 

with positive constants £, a and A . Here a will be a large 
even integer to accommodate a rolling potential having 
an asymptotic freedom for small 2 in the late Universe 



as GR with a cosmological constant Ao is recovered [21]. 
Current astronomical observation in relation to dark en- 
ergy suggests A = 10~ 122 ?~ 2 [22]. Since this value has 
negligible effects on inflation and core-collapse SNe, we 
shall neglect A for simplicity. 

From Eq. (2) we see that an asymptotic freedom of 
scalar coupling can arise from sufficiently large 2 or 
small T, for j3 < 0. Note T may also depend on im- 
plicitly, as demonstrated below. The dynamical balance 
between these two competing forces on the right hand 
side of Eq. (2) results in the trapping and relaxation 
of the scalar gravity in high energy environments (e.g. 
early Universe, core-collapse) and low energy environ- 
ments (e.g. in late Universe and solar system) regimes 
respectively. Apart from these guiding principles, the 
precise forms of 0(0) and U(0) are not critical. 

To focus on the salient effects, we shall neglect the 
spatial curvature parameter and associate the Einstein 
metric g a b with the Robertson- Walker spacetime with the 
Hubble parameter H(t) = a/a, the Einstein equation (1) 
yields the Friedmann equations 

H 2 = -r^r {U + u) 
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(U + 3T + U + 3p) 



with over-dot denoting time derivative. Here U and V 
are the effective energy density and pressure of de- 
rived from T ab in the standard manner and u = 4 u and 
p = £! 4 ]5 are the Einstein energy density and pressure 
of matter, in terms of the physical energy density u and 
pressure p of matter respectively, obtained from Eq. (3). 

As the equation of state for matter, we follow the stan- 
dard form p — w u for some constant w so that 

f= -u + 3p= -(l-3w)u. 

Since we do not expect matter to drive inflation, we shall 
be primarily interested in the range —1/3 < w < 1/3. 

After some algebra, the second Friedmann equation 
take the form: 

H + |(1 + w)H 2 + ^(1 - W )0 2 - c 2 (l + w)^4>" = 0. 

Under the trapping and quasistatic rolling of 0, we have 
« and w 0, and hence the scalar field equation (2) 
and the above equation become 

3(1 - 3w)H 2 p - 2c 2 (l - 3w)/3£0 CT + 2c 2 £cr0 ff ~ 2 = 

H + ^(1 + w)H 2 - c 2 (l + w)^ = 0. 
Furthermore, the trapped scalar satisfies 
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It is convenient to introduce Q(t) = q(t) + 1 using the 
decelerating parameter q(t), so that we can rewrite the 
above time evolution equation for the Hubble parameter 
in a neat form: 

using the following expression for Q(t): 

cr(l + w) 



Q(t) 



During inflation we have q(t) < and so Q(t) < 1, which 
through the above equation implies a approximate lower 
bound on (f> 2 : 

2 <t(3w + 1) 
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Once <fr 2 reaches (j) 2 . from above, inflation exits. The 
equivalent exit value for the Hubble parameter then fol- 
lows as 



H x 



There are two important limiting cases: As w — > 1/3, 
describing radiation, H x — > oo implying no inflation, con- 
sistent with cosmology where inflation preceding the ra- 
diation era. On the other hand, as w — > —1/3, violating 
the strong energy condition, H x — > and inflation never 
ends, as such exotic matter is capable of driving inflation. 

Suppose the scalar is initially trapped at time ti with 
value 0i = 4>(t{) with the corresponding Hubble parame- 
ter Hi = H(t\) > H x and Q[ = Q(ti), and then <j> slowly 
rolls towards zero. Then the time evolution of the Hubble 
parameter can be described by 



H (ty 



This in turn yields the approximate power law expansion 
for the scalae factor: 



a(t) wai [QiH^t - U) + 1] 



i/Qi 



Denote the exit time by i x so that H x = H(t x ) and a x = 
a(i x ). From the above relations we see that over the 
duration 

T = t x -t i *> Or 1 (#x 1 - Hr 1 ) 
the expansion of a(t) has the e-folding of 

lntox/aOwQ^Mfli/ifx). 

As a numerical example, if we take w = 0, j3 = 
-5, a = 20, i = 10- 22 ?- 2 Hi = 10-H- 1 then we see 
that the e-folding w 100 over the period r w 10 10 t p , 
which are compatible with the essential requirement for 
inflation. 



In the usual approach to inflation, it seems inconceiv- 
able to ask for observational effects, apart from imprints 
on cosmology such as CMB anisotropies and density per- 
turbations. However, since we are identifying inflation 
with scalar gravity, it is imperative to ask whether a re- 
lated scalar gradational effect would show up in astro- 
nomical events. As discussed above, such events neces- 
sarily involve high energy and density. Specifically, we 
shall address the core-collapse SN problem and assess 
how the ST theory relevant for inflation may count for 
(some of) the apparently missing energy required to ex- 
plain the observed powerful explosions [23]. 

To this end, we shall consider an additional contribu- 
tion to gravity due to in the Newtonian limit. We 
therefore approximate g a b by the Minkowski metric 77^ 
and assume a small scalar gravitational field. Therefore, 
we will ignore the potential V(<j>) as a 3> 1 and treat 4> as 
massless as per the asymptotic freedom of V(<j>) discussed 
above. We expect nontrivial effect to come from the cou- 
pling function fl(4>) in the presence of high density and 
adopt (3 < as in the inflation model above. 

To capture the essential dynamical trapping of <j> and 
its effect on gravitational binding energy, we consider an 
idealized spherical core-collapse model, in which the grav- 
itational collapse of an iron core of a heavy star is sud- 
denly halted as an incompressible neutron core of radius 
R with a constant (nuclear) density p is formed. From 
this moment, we shall establish the time evolution of the 
scalar field, with the associated energy shift and rele- 
vant timescale. We regard the collapse to be so rapid, 
it is not until the neutron core is just formed, that the 
scalar inside the core starts to evolve from a small cos- 
mological value towards a new equilibrium. However, the 
final equilibrium state of the scalar in our analysis below 
recovers features of "scalarization" discussed in the lit- 
erature [8-12]. We also note that the theoretical time 
dependent estimates below generate consistent numeri- 
cal simulation results in a different context reported in 
[24, 25]. Nonetheless, our analytical approach to the dy- 
namical trapping of scalar gravitational field offers fresh 
insight and simplicity for the present investigation. 

Given the high density of the neutron core, we shall 
neglect density outside the core. Using the ST theory 
in the Newtonian limit, the total gravitational potential 
depending on the radial distance r from the core center 
is derived from g 00 to be the sum <!> = $ n + $ s of the 
Newtonian gravitational potential $ n = (27r/3) Gp(r 2 — 
3R 2 ) and scalar gravitational potential $ s = (1/2) c 2 (3<j) 2 . 
Using the quasistatic pressure inside the neutron core 
satisfying Vp = — pV<I> and p = on the surface, we see 
that p(r) = — p[$(r) — &{R)]- By neglecting traveling 
scalar waves inside the neutron core, we represent <f)(r, t) 
for r < R as 



, s , x 2i? . 7rr 
<P(r,t) =<p{t)— sm — 



4 



where ip(t) = </>(0, t) denoting the value of the scalar field 
at the center of the neutron core. We can then substitute 
this into the scalar equation (2) using (3) and T = —c 2 p+ 
3p with p and p above, and average this equation over 
the volume of the neutron core. The approximate scalar 
energy inside the core can be calculated using T 00 to be 

E = l<fi 2 /2 + V(f) 

using the effective potential 

V(<p) =A<p 2 /2 + B<p 4 /4 

in terms of the following coefficients: 

I = (2/tt 2 )c 2 R 3 /G 

A = c 4 R/2G + ft (2.6 c 2 R 3 p- 7.7 GR 5 f) 
B = ft 2 (4.lc 2 R 3 p-8.9GpR 5 ). 

On the surface of the neutron core, the outgoing wave 
boundary condition applies. The power carried away by 
the scalar waves can be calculated using T° r to be 



4c 3 fl 2 
ir 2 G 



The conservation of scalar energy E + P = then yields 
the equation of motion for the amplitude <p as 



l(p + Aip + Bip 3 +C<p = 



(4) 



where C = (4/tt 2 ) c 3 R 2 /G. While B > is generally 
satisfied, A can be negative if ft is below the critical value: 



5.2 c 2 GR 2 p - 15 G 2 R 4 p 2 ' 

Then V(ip) becomes Higgs-like having local minima 
Vmin = -A/B with the minimum V(^ min ) = -A 2 IA.B. 
Since the relaxation time for Eq. (4) from ip w is — C / A 
the time to for ip(t) to evolve from a small cosmological 
value (f>o to (p m i n can be estimated using the relation 



A (po 

Let us now take again ft — — 5 as with the inflation 
example above, together with p = 3 x 10 17 kg/m 3 and 
R = 15 km. Then /3 C = —4.5 > ft and so scalarization will 
occur, leading to the dynamical trapping of scalar gravi- 
tational field in timescale to depending on the cosmolog- 
ical value 4>o- We see that to > 10 ms for <fi < 10~ 10 , a 
value compatible with current experimental tests. 

In this case, the trapped scalar has ipfain ~ ^-01 at 
the core center, and through V(</? m i n ), we find that this 
trapping releases 2.6 foe, or 0.15% of M Q c 2 , of gravita- 
tional energy, which is significant for a SN. This model 
therefore suggests that the discussed dynamical trapping 



and relaxation of scalar gravitational field could play an 
important role in explaining both inflation and delayed 
re-energization in core-collapse SNe. 

An opportunity potentially even exists to test the 
present model, since the amount of 44 Ti ejected from a 
SN, and thus available for detection by gamma-ray sen- 
sitive satellites [26], is thought to be a property of the 
underlying explosion mechanism. If the present mecha- 
nism were to affect the amount ejected a detection could 
discriminate between models. A recent experiment con- 
ducted at the CERN-ISOLDE facility [27] has attempted 
to remove the major nuclear physics uncertainty, enhanc- 
ing the chances of such a measurement being made. 
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